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ABSTRACT 

EK Draconis (HD 129333: G1.5 V) is a well-known yonng (50 Myr) solar 
analog. In 2012, Hubble Space Telescope retnrned to EK Dra to follow np a 
far-nltraviolet (FUV) SNAPshot visit by Cosmic Origins Spectrograph (COS) 
two years earlier. The brief SNAP pointing had fonnd snrprisingly redshifted, 
impnlsively variable snbcoronal “hot-line” emission of Si IV 1400 A (T ~ 8 x 
10^ K). Serendipitously, the 2012 follow-on program witnessed one of the largest 
FUV flares ever recorded on a snnlike star, which again displayed strong redshifts 
(downflows) of 30-40 km s“^, even after compensating for small systematics in the 
COS velocity scales, uncovered through a cross-calibration by Space Telescope 
Imaging Spectrograph (STIS). The (now reduced, but still substantial) ~10 km 
s“^ hot-line redshifts outside the flaring interval did not vary with rotational 
phase, so cannot be caused by “Doppler Imaging” (bright surface patches near a 
receding limb). Density diagnostic O IV] 1400 A multiplet line ratios of EK Dra 
suggest He ~ 10^^ cm“^, an order of magnitude larger than in low-activity solar 
twin a Centauri A, but typical of densities inferred in large stellar soft X-ray 
events. The self-similar FUV hot-line prohles between the flare decay and the 
subsequent more quiet periods, and the unchanging but high densities, reinforce a 
long-standing idea that the coronae of hyperactive dwarfs are flaring all the time, 
in a scale-free way; a flare-ona if you will. In this picture, the subsonic hot-line 
downflows probably are a byproduct of the post-flare cooling process, something 
like “coronal rain” on the Sun. All in all, the new STIS/COS program documents 
a complex, energetic, dynamic outer atmosphere of the young sunlike star. 

Subject headings: ultraviolet: stars — stars: individual (HD 129333= EK Dra; 
HD 128620= a Cen A) — stars: coronae 
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1. INTRODUCTION 

Flares on the Sun result from catastrophic releases of energy previously stored in strong 
tangled surface magnetic helds. These outbursts are most conspicuous at high energies, and 
under the right circumstances (or wrong, depending on your perspective) can profoundly 
impact the heliosphere and Earth. This is one reason why flares have drawn intense interest 
from solar observers and theorists for many decades; and indeed have spawned a vibrant 
Space Weather Community in recent years. A less appreciated aspect of solar flares involves 
the stellar cousins, where extraordinarily large X-ray (and even white-light) outbursts have 
been recorded, in some cases dwarfing by a considerable margin the largest known flares 
on our own Sun, albeit a middle-aged and magnetically rather uninspired star. The most 
extreme of these are “superflare^,” which at high energies can exceed the quiescent coronal 
X-ray emission of a star by factors of a hundred or more. Interest in such rare, but poten¬ 
tially catastrophic, events concerns how they might affect, for example, survival of primitive 
lifeforms struggling for existence on a semi-habitable world orbiting a young, magnetically 
active star. 

Until recently, reports of superflares mainly involved interpretations of historical ma¬ 
terials, such as long-term astronomical plate collections, and disconnected isolated reports 
or inferences (Schaefer et al. 2000). Now, however, the Kepler Mission has provided more 
substantive statistics on superflares - although limited to white-light counterparts - from 
a wide variety of stars, normal and otherwise (e.g., Maehara et al. 2012; Shibayama et al. 
2013; Candelaresi et al. 2014). Further, the Swift gamma-ray burst satellite has witnessed 
a handful of the most extreme events, radiating in hard X-rays, thanks to its broad view 
of the high-energy sky and rapid pointing response (e.g., Osten et al. 2010; Drake et al. 
2014). The Swift alerts also have promoted supporting multi-spectral observations of flares, 
especially at radio frequencies that capture non-thermal emissions from accelerated particles 
(e.g.. Fender et al. 2015). On the other hand, there have been few cases of such outbursts 
recorded at wavelengths where the powerful diagnostics of high-resolution spectroscopy can 
be brought to bear, namely the ultraviolet. The lack of good examples mostly is because 
these rare events are, well, rare. 

Serendipitously, however, a recent HST Cosmic Origins Spectrograph (COS) campaign 
on young (~50 Myr) solar analog EK Draconis (HD 129333; G1.5 V), well-described in the 
previous literature (e.g., Konig et al. 2005, and references therein), caught a large hour-long 
FUV outburst, in “hot lines” like the C IV 1550 A doublet (T~ 1 x 10® K) and the extremely 
hot [Fe XXI] 1354 A coronal forbidden line (~ 10^ K). The EK Dra flare, viewed with 


^Originally called “flashes” by Schaefer (1989), who first called attention to the class. 
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the excellent spectral resolution, broad wavelength coverage, and high sensitivity of COS, 
provides a unique opportunity to explore properties of such events, to inform interpretations 
of the physical processes involved. Moreover, the behavior of the FUV spectrum outside 
the flare, over the several days of intermittent attention by HST, opens a window onto the 
“quiescent” properties of a hyperactive chromosphere and corona, again with unprecedented 
clarity. 


2. OBSERVATIONS 
2.1. Previous HST Studies 

There were two previous HST spectroscopic programs of note on EK Dra. In 1996, Saar 
& Bookbinder (1998) carried out high time resolution measurements with the low-resolution 
FUV grating of Goddard High-Resolution Spectrograph (GHRS) over a continuous interval 
of about six hours (about 10% of the spin period). The authors identihed several factor of 
2-3 impulsive flare-like bursts in Si IV 1400 A and C IV 1550 A when the light curves were 
binned into 30 s intervals. The identihable hot-line emission enhancements occupied about 
8% of the timeline. The remainder was attributed to “quiescent” periods. 

A decade and a half later, in 2010, Ayres & France (2010) targeted EK Dra with newly 
commissioned COS as part of a larger “SNAPshot” (partial-orbit Eller observations) survey 
of FUV coronal forbidden lines, such as [Fe XXI] A1354 (T ~ 10^ K), in more than a dozen 
late-type stars. The objective was to measure Doppler profiles of high-energy coronal lines at 
FUV resolution, some 20 times better than currently possible in the soft X-ray band, the more 
natural - but spectroscopically challenged ~ home of coronal transitions. The COS SNAP 
observation of EK Dra covered only a short time span, about 20 minutes, compared with the 
earlier GHRS pointing. But, it had the important advantage of achieving moderately high 
spectral resolution (~ 20 km s“^) in the integrated spectrum, as well as comparably high 
time resolution (30 s) for flux measurements of key species like C II 1335 A (2x10^ K), Si IV, 
and [Fe XXI] collected in the single G130M grating setting feasible for the brief partial-orbit 
visit. 


The Si IV doublet, in particular, exhibited twin factor of ~2 enhancements over the 
mere 20 minute interval, remarkably without comparable responses at lower temperatures 
(C II) or higher (Fe XXI), contrary to expectations for normal flares, which tend to impact 
all the atmospheric temperature layers more-or-less simultaneously. Equally remarkable, the 
Si IV doublet features displayed quite large redshifts, about 20 km s“\ relative to narrow 
chromospheric features (which set the empirical zero-point velocity for the COS exposures). 
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Redshifts of the hot lines are well documented on the Sun, attributed to nearly radial 
downflows of subcoronal material, but the average shifts in a disk-integrated spectrum might 
be only a few km s“^ or less (Peter 2006). Ayres & France interpreted the large redshifts of 
Si IV on EK Dra as a signature of a scaled up version of solar “coronal rain;” previously hot 
multi-MK gas suspended high in magnetic loops rapidly cooling and draining back down to 
the lower atmosphere, perhaps ~ in the case of EK Dra - in a catastrophic way (what one 
might call an anti-flare: impulsive cooling rather than heating). 


2.2. The 2012 HST Campaign 

The follow-up HST program, described here, partnered the high sensitivity of COS 
with the excellent wavelength precision of STIS to explore alternative explanations for the 
COS redshifts, and thus also the coronal rain hypothesis: (1) subtle distortions in the COS 
wavelength scales; and (2) “Doppler Imaging” effects, whereby active patches on the stellar 
surface might appear in the integrated spectrum as a redshift if, say, a bright region happened 
to be near the receding limb (where the projected rotational velocity could be as much 
as -|-nsini ~ -\-17 km s“^). A secondary objective was to carry out a more extensive 
time-domain study of FUV variability, along the lines of the Saar & Bookbinder work, but 
with the added dimension of time-resolved spectroscopy possible with high-sensitivity, time- 
tagged COS. A related goal was to extend the previous minimal SNAP wavelength coverage 
(G130M/detector side A: 1290-1430 A), adding G160M/sides A and B, especially to reach 
the important C IV 1550 A doublet and He II A1640 at the longer wavelengths; and with 
STIS to probe below the COS G130M/A shortwavelength cutoff (side B was deactivated 
because H I Lya A1215 violated detector bright limits), including especially Lya itself. 

The characteristics and operational capabilities of COS have been described in a number 
of previous publications, especially Green et al. (2012); and companion UV spectrograph 
STIS, by Woodgate et al. (1998) and Kimble et al. (1998). 

The new HST program was carried out over two days in late-March 2012 (the stellar 
rotation period is 2.8 d). EK Dra was in the HST Continuous Viewing Zone (CVZ) at the 
time, by design, which allowed each of the six allocated 96 minute orbits to be used fully, 
without any loss to Earth occultations. Figure 1 is a schematic timeline of the program, also 
summarized in Table 1. 
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Fig. 1.— Observational timeline of the 2012 EK Dra HST program. Circled numerals are visit 
numbers. Yellow blocks represent initial STIS E140M exposure, split into four equal length sub¬ 
exposures; preceded by a target acquisition and peak-up (thin blue vertical lines). Eirst COS 
observation (visit 2) followed immediately, with grating setting G160M, then G130M, both split 
into four EP-POS steps. Each EP-POS sub-exposure is represented by a streak image derived from 
the associated event list (darker shading represents higher count rates). Non-uniform spacing of 
the G130M sub-exposures in visits 2 and 5 devolve from the scheduling process. Vertical red lines 
mark the COS NUV imaging target acquisitions. 
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Table 1. HST Observations of EK Dra (2010-2014) 


Dataset 

UT Start 

^exp (s) 

Aperture 

Grating-Acen(A) 

FP-POS 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


2010 SNAP Program 

Ib3e34010 

Apr-22 09:18 

1160 

PSA 

G130M-1291 

3,4 

2012 Visits 1 & 2 (contiguous) 

oboqOlOlO 

Mar-27 16:41 

8880 

PHT 

E140M-1425 


lboq02010 

Mar-27 19:57 

1000 

PSA 

G160M-1577 

3,4,1,2 

lboq02020 

Mar-27 20:41 

1780 

PSA 

G130M-1291 

3,4,1,2 

2012 Visit 3 

IboqOSOlO 

Mar-28 05:43 

1000 

PSA 

G160M-1577 

3,4,1,2 

lboq03020 

Mar-28 06:48 

1780 

PSA 

G130M-1291 

3,4,1,2 

2012 Visit 4 

lboq04010 

Mar-28 20:33 

1000 

PSA 

G160M-1577 

3,4,1,2 

lboq04020 

Mar-28 21:16 

1780 

PSA 

G130M-1291 

3,4,1,2 

2012 Visit 5 

lboq05010 

Mar-29 10:25 

1000 

PSA 

G160M-1577 

3,4,1,2 

lboq05020 

Mar-29 11:08 

1780 

PSA 

G130M-1291 

3,4,1,2 

2014 SNAP Program 

ocd839010 

May-06 10:10 

600 

PHT 

E230H-2713 


Note. — 

Col. 1 prefix “o' 

” is STIS, 

“1” is COS 

1. Gol. 4 “PHT” 

is the STIS 


0.2" X0.2" photometric aperture; PSA is the COS 2.5"-diameter Primary Science 
Aperture. Col. 6 indicates sequence of (standard) small grating steps to mitigate 
fixed pattern noise. The 2014 STIS observation was from an unrelated program 
(S. Redfield P.I.), but nonetheless relevant to the present study. 

























7 


The first visit, consisting of two orbits, featnred a deep STIS exposure with the medium- 
resolution FUV echelle (E140M-1425) for 8.8 ks, split into four equal-length sub-exposures (to 
mitigate instrumental drifts due to telescope “breathing”). The initial target acquisition was 
with the CCD in visible light through the F25ND3 hlter. Following that, a peak-up (accurate 
target centering via a raster search) was performed with the 0.2"x0.06" slit in dispersed 
visible light (CCD with G430M-3936 grating setting). Then, the STIS FUV exposures were 
taken, through the 0.2"x0.2" photometric aperture for maximum throughput, although in 
principle retaining high accuracy in the wavelength scales thanks to the narrow-slit peak-up. 

Immediately following the STIS segment, although in a separate visitEI, was the hrst of 
four one-orbit COS pointings, spaced evenly over roughly half a spin cycle so all sides of the 
possibly inhomogeneous surface activity would be visible. 

The four COS visits ostensibly were identical and occupied a single CVZ orbit each. The 
COS visits began with an NUV imaging target acquisition, using Mirror-A and the Bright 
Object Aperture (BOA). Following that, a 1000 s observation was taken with G160M-1577 
(1387-1557 A on detector side B; 1578-1749 A on side A), divided into equal length sub¬ 
exposures, one at each of the standard four “FP-POS” positions (small grating rotations to 
suppress hxed-pattern noise). Then, a 1780 s integration, also split into four equal FP-POS 
sub-exposures, was taken with G130M-1291 (1292-1432 A on side A; side B turned off, as in 
the earlier SNAP program, owing to overbright Lya). 

All the COS observations were done in TIME-TAG mode, which delivers a time-stamped 
record of photon (x, y) positions in the 2-D detector format (x corresponds to wavelength; y 
to cross-dispersion spatial coordinate). These photon event lists can be processed in a number 
of ways to provide time histories of fluxes and/or line prohles. The particular G130M and 
G160M settings were chosen because they overlap at the important Si IV 1400 A doublet, 
so that it would be recorded in all the observations. Unfortunately, the key [Fe XXI] coronal 
forbidden line was captured only in the G130M segments; and the important subcoronal 
G IV doublet only in the G160Ms. 


2.2.1. STIS Observations in Visit 1 

The initial STIS observation and hrst of the GOS visits were intended to be carried out 
as a pair, as close together as practical, to serve as a velocity cross-validation. Although GOS 


^COS utilizes different Fine Guidance Sensor (FGS) sky segments (“pickles”), and thus a different set of 
Guide Stars; so a separate Guide-Star/target acquisition had to be performed. 



exposures have an internal wavelength reference (lamp flash), the calibration spectrum falls 
well above the stellar stripe, on a part of the detector where uncompensated geometrical 
distortions are different than where the stellar spectrum is located. Further, the specihc 
^/-position of the stellar spectrum, and thus the local geometrical distortions it experiences, 
depends on the acquisition method. However, because the same strategy was used for all 
four COS visits, any wavelength issues identihed in the comparison of the initial COS visit 
with STIS should apply equally well to the three subsequent COS pointings. At the same 
time, while STIS has exquisite velocity precision and can record S/N~20 prohles of the 
brightest lines (including Lya, forbidden to COS), it lacks the sensitivity to resolve shapes 
of the bright lines on short time scales, or dig down to fainter features like the key O IV] 
1400 A density diagnostic (see, e.g.. Cook et al. 1995, and references therein). 

The STIS spectra of EK Dra were processed and combined using protocols developed 
for the HST Advanced Spectral Library Project (ASTRAlj^), which includes a post-pipeline 
correction for subtle wavelength distortions. The four independent sub-exposures were co¬ 
added, but without the usual cross-correlation registration, because the relatively faint spec¬ 
tra lacked narrow features of suitable strength for the purpose. Figure 2 compares represen¬ 
tative intervals of the average STIS FUV echelle spectrum of EK Dra to low-activity solar 
twin a Centauri A (G2 V), also observed by STIS with the same setup (see Ayres 2015). 

At hrst glance, the two spectra are surprisingly alike, despite the enormous gulf between 
the stars in coronal activity (the Lx of EK Dra is about a thousand times that of a Cen). For 
example, the narrow chromospheric lines of O I, C I, and Cl I are quite similar in appearance 
between the hyperactive star and low-activity counterpart. Even the hot lines like C IV have 
similar emission cores. 

Closer examination, however, reveals several differences. The most noticeable is the 
exceptionally broad H I Lya prohle of EK Dra (although somewhat emphasized by the 
normalization approach: more distant EK Dra [d ~ 34 pc] is more heavily absorbed by 
interstellar H I than close-by a Cen [1.3 pc], so is missing much of its intrinsic narrow peak). 
A second, in fact blatant, difference is the bright [Fe XXI] coronal forbidden line of EK Dra, 
which is absent in low-activity a Cen (the weak, narrow feature at 1354 A is a chromospheric 
C I line). On the other hand, the striking difference at the He II A1640 Balmer-a line mainly 
comes about because low-activity a Cen is dominated by a chromospheric Fe II blend on the 
blueward side of He II; rather than high-excitation He II itself as in EK Dra. The prominent 
photospheric continuum of a Cen A in the He II A1640 panel is an artifact of the scaling 
procedure: the line-to-continuum contrast is much smaller in the low-activity star compared 


^see: http://casa.Colorado.edu/^ayres/ASTRAL/ 
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Fig. 2.— Examples of STIS E140M FUV echelle spectra of EK Dra (open circles; 4-point smoothing 
[2 resolution elements (resels)]) compared to low-activity solar analog a Cen A (gray shading 
outlined by black dots). The x-axis scale in each panel is A A (A) relative to the wavelength cited 
in the axis title. Red dotted-dashed curves are the smoothed 1 a photometric noise levels (per 
resel) for EK Dra. The two spectra were adjusted to match the peak of one of the spectral features 
in each interval. The wavelength scales were registered to the average velocity shift of narrow 
chromospheric emissions. Conspicuous absorptions in several of the chromospheric features (e.g., 
H I A1215 Lya, O I A1302, and C II A1334) are from the interstellar medium (ISM). The ISM 
velocity is more redshifted in EK Dra, relative to the stellar radial velocity, than in a Cen A, so 
the EK Dra absorptions cut into the red edges of the narrower emission cores (e.g., O I A1302 and 
C II A1334). The small peak in the center of the EK Dra Lya profile is atomic hydrogen skyglow 
through the 0.2" x 0.2" aperture. The dip in the blue wing of Si III A1206 in a Cen is an artifact due 
to contamination of the local background by Lya light scattered from the adjacent echelle order. 
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with hyperactive EK Dra. If the normalization instead had been according to the bolometric 
fluxes, the 1640 A photospheric continuum levels of the two stars would be similar, but the 
He II emission of EK Dra would dwarf that of a Cen A. (Indeed, a similar dwarhng would 
occur in all the panels with an /boi normalization.) More subtle, but significant, differences 
are the broad low-intensity wings of the EK Dra hot lines: Si III, Si IV, C IV, and N V 
(T ~ 0.5 — 2 X 10^ K). There are some additional differences seen in the 1353.7 A panel, 
having to do with pumping of the Cl I and C I lines by higher excitation species, leaving O I 
A1355 much weaker by comparison in the active star. 


+ 100 
+50 
0 

-50 


Fig. 3.— Schematic EK Dra COS FUV spectrum from the 2012 campaign. The G130M and 
G160M settings overlap in the important 1400 A region. Note the detector gap between G160M 
sides A and B at 1560-1580 A. Image was constructed from the totality of the associated photon 
event lists from GOS visits 2-5, with adjustments in the detector y centroids to align the three 
separate segments, for display purposes. 



2.2.2. COS Observations in Visits 2-5 

Figure 3 is a schematic of the COS FUV spectrum covered by the 2012 program. The 
lamp-flash calibrations (barely visible at this intensity stretch) can be seen 100 pixels above 
the main stellar stripe. The two G160M segments were observed simultaneously. The long- 
wavelength end of G130M on detector side A overlaps the shortwavelength portion of G160M 
on side B at the key Si IV 1400 A resonance doublet (and O IV] density-diagnostic multi- 
plet); side B of G130M was disallowed, as in the 2010 SNAP program, owing to intense 
chromospheric Lya. The broad features (in cross-dispersion y coordinate) at 1300 A are the 
O I resonance triplet, dominated by sky-glow through the 2.5"-diameter GOS aperture. The 
other bright features are G II A1335, Si IV A1400 (distinct pair), and G IV A1550 (close pair). 














11 


Table 2. Measurement Parameters for the COS Event Lists 



Integrated Count Rates 



Spectral Profiles 

Type 

Band (A) 

(S/N)erit 

Ao (A) 

AA (A) 

Acrit (counts) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


G130M-1291 

LINE-1 

1335.1±1.5 

20 

1335.703 

0.05 

3000 

LINE-2 

1354.0±0.75 

10 

1354.080 

0.10 

1000 

LINE-3 

1393.8±1.0±1402.8±0.75 

20 

1402.773 

0.05 

1000 

GON-A 

1344.5±6.0±1381.0±8.0 

20 




BKG-A 

1360±70 





G160M-1577 

LINE-1 

1393.8±1.0±1402.8±0.75 

20 

1402.773 

0.05 

1000 

LINE-2 

1549.5±2.75 

20 

1548.204 

0.05 

2000 

LINE-3 

1640.35±1.0 

20 

1640.404 

0.05 

1000 

GON-B 

1435±25 

20 




BKG-B 

1475±85 





GON-A 

1610±25 

20 




BKG-A 

1665±85 






Note. — Col. 1: “LINE” for emission line, “CON” for continuum band, “BKG” for 
off-spectrum background. Suffixes “B” and “A” refer to the respective detector segments. 
For the integrated count rates, col 2. is the integration bandpass, and col. 3 is the signal- 
to-noise of each measurement (defining the duration. At, of the bin). For the line shape 
extractions, col. 4 is the central wavelength of the constructed prohle (±250 km s“^ 
bandpass), col. 5 is the bin size, and col. 6 is the number of counts collected in each profile 
(again, defining a At). For the fluxes and profiles, the cross-dispersion (j/) extraction 
window was ±9 pixels. Background swaths were positioned at Ay = ±40 pixels from the 
spectrum center, with 20 pixel widths. Background count rates were determined from the 
full duration of each subexposure, regardless of S/N. 










As mentioned earlier, the initial STIS and COS visits were intended to be as close as 
possible to enable a cross-comparison of the two UV instruments. Thanks to diligent work 
by the HST schedulers, the first COS exposure (G160M) began a mere 40 minutes after 
the hnal STIS readout, normally close enough to ensure comparable stellar conditions over 
the three orbits of contiguous STIS and COS attention. Indeed, the STIS FUV emission 
line intensities were stable over the 2.5 hours of near-continuous E140M exposure. As luck 
would have it, however, less than an hour later at the beginning of the first COS G160M 
sequence, EK Dra apparently was in outburst, experiencing a giant FUV flare. The star had 
brightened up in C IV by an order of magnitude, nearly triggering the COS detector bright 
limits, a remarkable performance for a nearly 8th-magnitude G dwarf. 

Figure 4a illustrates time-resolved line fluxes and prohles from the hrst COS observation 
(visit 2), with grating G160M. The top panel displays the total spectral image (all four FP- 
POS event lists combined), with specihc integration bandpasses marked for the Si IV doublet 
(blue), C IV (red), and He II (green); as well as two line-free continuum bands (“B:” orange 
dashed; “A:” orange/black dashed); and two pairs of flanking background swaths (black 
dashed, and black-yellow dashed), one for each detector segment. The parameters dehning 
these measurements, and the subsequently described spectral prohle extractions, are listed 
in Table 2. 

The middle panel of Fig. 4a compares the time-resolved line and continuum count rates, 
with the same color-coding. The individual points represent counts accumulated in sliding 
time bins to achieve a specihed signal-to-noise ratio (S/N)crit: see Table 2. The two contin¬ 
uum bands (orange, and black outlined orange, circles) were scaled down a factor of 4 for 
display purposes, while the background count rates (black and yellow dots; taken from each 
full FP-POS time interval) were scaled upward by 100 times relative to the continuum bands 
(taking into account the continuum integration wavelengths and cross-dispersion width; see 
Table 2) to illustrate the possible influence of the background for the worst-case scenario 
(broad wavelength bands). The background corrections for the continuum bands in this 
instance are negligible, but still were compensated. The background contributions for the 
tiny line spectral footprints are entirely negligible, and were not corrected. 

The lower panels illustrate spectrally resolved prohles of three bright features - Si IV 
A1402, He II A1640, and C IV A1548 - each accumulated in time intervals corresponding 
to a hxed number of counts {Ncrit- see Table 2). Si IV A1402 was chosen over stronger 
Si IV A1393, because the latter is affected by bad pixels in one of the FP-POS segments of 
G160M/B. The specihc accumulation intervals are highlighted by small ticks in the middle 
panel time sequences. In each of the bottom spectral panels, the time-resolved prohles are 
the thin colored curves. The background yellow shaded envelope represents a prohle of the 
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same transition from a spectrum averaged over visits 3-5 (see below), scaled to the maximum 
and minimum of the collection of line peaks for the particular visit. The upper and lower 
outlines of the envelopes are black dotted. Remarkably, despite the large flux changes during 
the post-flare decline caught by visit 2, the prohles of all three features are nearly identical to 
the scaled averages; showing the characteristic redshift of the line core, and the high-velocity 
extended wings, especially on the redward side of line center. 

Figure 4b is similar, for the G130M side A spectrum in the second part of visit 2. Here, 
the integrated count rates and extracted prohles are for C II, [Fe XXI], and Si IV. The two 
continuum bands (now smaller, to avoid the more numerous chromospheric emissions at 
these wavelengths) were combined into a single continuum hux, scaled down a factor of two 
in the hgure for display. The background count rates, again adjusted to the total (Ax, Ay) 
continuum bandpass, were scaled upward a factor of 25 relative to the displayed continuum 
values. All the huxes continue to decline during this interval, but again the time-resolved 
line prohles are self-similar to the visits 3-5 averages. 

Figure 5a illustrates the G160M sequence of visit 3, about 10 hours after the initial hare 
segment. The continuum and background count rates are scaled as in Fig. 4a. All the huxes, 
especially the continuum bands, have fallen compared with visit 2, and now the background 
corrections are more signihcant. The line prohles span a more compressed range, essentially 
identical to the visits 3-5 average. 

Figure 5b depicts the G130M exposures of visit 3. The continuum and background 
huxes are scaled as in Fig. 4b. Again, there is little variability in the integrated huxes or 
the time-resolved prohles. The analogous diagrams for later visits 4 and 5 are very similar 
to those of visit 3, and are not illustrated. 
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Fig. 4a.— Time-resolved line fluxes and profiles from the first COS observation (visit 2), during 
the large FUV flare, with grating G160M (see text for detailed description). 



Fig. 4b. 


Same as Fig. 4a, for the G130M side-A exposures during second part of visit 2. 

































































Figure 6 is an expanded view of the flare in Si IV, which was recorded in both COS grat¬ 
ing settings. As in Figs. 4 and 5, the symbols represent sliding time bins corresponding to 
(S/N)crit = 20. Accordingly, the bins become wider in time, and the points more separated, 
as the count rates decline during the later stages of the event. The COS Si IV values were 
normalized to the average of visits 3-5, which were much quieter than the visit 2 flare. The 
G160M/B and G130M/A Si IV time series were separately normalized to the corresponding 
quiescent averages, to avoid possible differences in sensitivity with the different grating set¬ 
tings on the different detector segments. The STIS Si IV points from visit 1, immediately 
prior to the flare rise, were based on calibrated fluxes integrated over the doublet compo¬ 
nents in each of the four independent E140M sub-exposures, divided by the integrated flux, 
measured in the same way, from an average COS spectrum constructed from the quiescent 
visits 3-5 (as detailed later). The flare evolution was modeled by twin exponentials: a rapid 
initial decline with an e-folding time constant of 7 minutes, on top of a slower decay with 
a 1.9 hour time scale. Such bimodal decays are common among stellar X-ray flares, for 
reasons that are not completely understood (see examples and discussion by Tsang et al. 
2012). The lower panel of Fig. 6 provides a broader view of the Si IV behavior during the 
entire HST program, showing explicitly the relatively quiet intervals prior to and following 
the flare event. 

Figure 7 is a multi-spectral view of the EK Dra flare decay, and subsequent quiescent 
intervals. Error bars to the left of the vertical dotted line indicate the STIS average levels 
and standard deviations (over the four sub-exposures). The COS time series to the right are 
the same as illustrated in Figs. 4 and 5, but now including visits 4 and 5, and normalized 
to the visits 3-5 average. Note the small flare decay that affects all the features in the 
second part of visit 3, together with a number of Si IV and C IV transients like those seen 
in the Saar & Bookbinder time series and the 2010 SNAP program. Also note the two 
bursts (at t =3.6 hr and 3.7 hr) in the G160M blue and red continuum bands (marked by 
the corresponding colors), which remarkably are not repeated in any of the other activity 
tracers. This phenomenon is discussed in more detail later (§3.4). 

Although the large FUV flare seemingly ruined the carefully orchestrated STIS-COS 
velocity cross-calibration, the bright side is that this was the first such signihcant event on 
a sunlike star to have been captured with the un-matched spectroscopic and time-resolution 
capabilities of an instrument like COS; a unique opportunity to explore properties of the 
outburst beyond the usual simple intensity evolution (with modest temperature discrimi¬ 
nation) delivered by, say, an X-ray telescope. And, as described above, following the large 
flare EK Dra settled into an period of apparent quiescence - to the extent that a star with 
a thousand times the X-ray luminosity of the Sun can be considered “quiet” - for the next 
three COS visits (spanning a day and a half, corresponding to about half a rotation period). 
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Fig. 5a.— Same as Fig. 4a for the G160M sequence of visit 3, about 10 hours after the initial 
flare segments. Unlike Fig. 4a, there is little variability of the integrated fluxes or the time-resolved 
profiles. 
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Fig. 5b. 


Same as Fig. 4b for the G130M sequence of visit 3. 
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Fig. 6.— Si IV A1393 + A1402 time series covering visits 1 and 2 (top panel), and the full campaign 
(bottom panel). Origin of the timeline is the initial STIS target acquisition. The Si IV count rates 
were normalized to the average of visits 3-5. Right side of the upper panel contains a close-up view 
of the serendipitous large FUV flare, captured in Si IV by the sequential grating settings (each split 
into four FP-POS sub-exposures). The flare decay can be described as an initial rapid exponential 
drop on the scale of minutes (green dotted-dashed curve), followed by a slower exponential decline 
on the scale of hours (blue dotted-dashed; red dashed is sum). Left side of the upper panel depicts 
the four STIS sub-exposures of visit 1. Lower panel portrays the 1.5-hour event within the context 
of the pre-flare (STIS visit 1) and post-flare (COS visits 3-5) quiescent periods. 
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Fig. 7.— Multi-spectral view of the EK Dra flare decay, and subsequent quiescent intervals. Error 
bars to left of the dotted line (which marks beginning of visit 2) indicate the STIS average levels. 
Visit numbers are in red in the upper panels. Count rates of each species were normalized to the 
visits 3-5 average. Larger gray symbols in each panel (except top-most) are the Si IV time series 
repeated for comparison. 
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The similarity of the time-averaged FUV spectra of these three visits, both among them¬ 
selves and compared with the initial STIS pointing, ultimately provided the desired velocity 
cross-validation (as described later in §3.1). So, all was not lost. 

Furthermore, the lack of significant FUV prohle variability of the hot lines between 
STIS visit 1 and COS visits 3-5, implies that “Doppler Imaging” was not an important 
factor during the 2012 HST campaign. In fact, the constancy of the hot-line fluxes over 
the more than half a rotational cycle also indicates that the surface distribution of FUV- 
emitting regions was relatively uniform, at least in that epoch. At higher energies, in contrast, 
Gudel et al. (1995) reported factor of ~2 systematic variations in ROSAT all-sky survey X- 
ray fluxes, in step with the rotation period, and a suggestion of correlated variability at 
microwave frequencies as well. On the Sun, albeit perhaps not the best guide to hyperactive 
EK Dra, the FUV-bright plage regions surrounding sunspot groups are much more spatially 
pervasive than the dark spots themselves, and more so than the compact coronal magnetic 
loop systems associated with the active region as well. 


2.2.3. EK Draconis B 

EK Dra is reported to have a faint, M-type companion in a highly eccentric 45 year 
orbit (Konig et al. 2005). Based on the authors’ ephemeris, the secondary would have been 
just past maximum elongation in 2012, nearly due south of the primary with a separation of 
0.75". This object would be at the inner edge of the COS PSA, and potentially visible as a 
slightly displaced spectrum (unfortunately on top of the main spectrum in the +x direction, 
given the position angle of the COS focal plane at the time of the 2012 program), but entirely 
outside the STIS 0.2"x0.2" slot. However, there is no obvious line doubling due to a second 
spectrum visible in any of the COS FUV grating images; although the predicted separation 
between the two is only -1-0.35 A (EK Dra B to red of A), and the bright primary lines are 
broad. 

The secondary also was not clearly visible in (stacked versions) of either the STIS or 
COS acquisition images of the 2012 campaign. Serendipitously, however, the faint companion 
was recovered in a subsequent STIS pointing on EK Dra, in an 2014 ISM SNAP program (S. 
Redheld P.I.; see Table 1). The partial orbit visit had a more conservative CCD acquisition: a 
total exposure time of 2.3 s through the F25ND3 aperture, signihcantly deeper than the 0.3 s 
of the 2012 program with the same hlter. The secondary was about 1% of the primary’s 
brightness, at a position angle of 172° and separation of 0.75", just as predicted by the 
ephemeris of Konig et al. (2005). The V magnitude of the secondary would be about -1-13.5, 
taking into account that the broad-band CCD response favors red stars over yellow stars. 
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The primary/secondary AV ~ 6 is what was suggested by Konig et ah 

This is important, because a youthful red dwarf could be quite active, and its hot 
lines, if bright enough, potentially could affect the prohles of the primary features (with the 
unfortunate location of B on top of the A spectrum in 2012). However, scaling FUV fluxes 
of the most active dMe flare stars, like AU Microscopii and AD Leonis, to the estimated V 
magnitude of the secondary predicts peak flux densities, say at C IV, of only a few percent of 
the primary’s, which would be of little concern. A negligible impact by the red star, at least 
on the quiescent FUV spectrum of EK Dra, is supported by the fact that the COS visits 3-5 
spectra are so similar to STIS visit 1, which entirely excluded the secondary by virtue of the 
small STIS aperture. 

As for the large FUV flare, if it were from the red star, the associated soft X-ray 
luminosity (see below) would have been quasi-bolometric (i.e., Lx/Tboi ~ 1); not impossible, 
but very unlikely. The ~2 known examples, such as the quasi-bolometric flare on DG CVn 
reported by Drake et al. (2014), have come from the Swift Burst Alert Telescope, and soft- 
X-ray follow-up by that satellite. The survey burst alert situation is much more likely to 
capture extremely rare events because it is viewing a whole sky full of objects, all the time, 
rather than pointed at a single target for a limited period. 

Furthermore, the enhanced narrow peaks of the FUV hot lines during the flare were 
nearly identical in shape and velocity to the quiescent prohles of the G star seen outside the 
outburst; rather than shifted by the ~-|-75 km s“^ to the red expected if the FUV lines were 
instead from the secondary star. The preponderance of evidence is that the large FUV hare 
was from the G-type primary. 


3. ANALYSIS 

3.1. STIS-COS Wavelength Cross-Calibration 

A main objective of the 2012 EK Dra campaign was to test the reality of the large 
~ 20 km s“^ Doppler shifts of Si IV recorded in the brief COS SNAPshot two years earlier. 
The core of that test was to validate the COS wavelengths against the highly precise scales 
of STIS. In particular, the E140M-1425 echelle setting captures the entire EUV region (1150- 
1700 A) in a single shot, so there is no issue trying to tie together multiple spectral segments 
(e.g., G130M and G160M of COS, including the individual FP-POS steps, each with its own 
potential calibration issues). Further, the STIS MAMA cameras have hard-wired physical 
pixels, whereas the cross-delay-line detectors of COS have floating pixels dehned by the 
{x, y) impact point of the electron cloud associated with a photon event. The implication for 



STIS is that a given pixel always experiences whatever geometrical distortions are present 
locally, and these can be accounted to some degree in the wavelength dispersion solution. 
With COS, conversely, the geometrical distortions are sampled quasi-randomly by whatever 
distribution of impact points happen to fall at a particular location. 

Perhaps more important, the STIS wavelength calibration light (from onboard hollow- 
cathode emission lamps) follows the same optical path as an external stellar point source, 
falling on the same echelle order pixels as the stellar spectrum, and consequently experiences 
the same geometrical distortions. Thus, the STIS wavelength calibration spectra encode 
signihcant information concerning these local geometrical distortions, which then can be 
accommodated in a hybrid dispersion solution, such as described in the ASTRAL protocols 
mentioned earlier. For COS, on the other hand, the lamp-flash spectrum falls well above 
the spectral stripe (see, e.g.. Fig. 3), on a portion of the detector where the geometrical 
distortions are completely different. Consequently, the COS lamp-flashes are less relevant 
for conhrming the pipeline dispersion solutions. 

To test the wavelength hdelity of COS, the three independent combinations of grating 
setting and detector segment were considered separately. Reference spectra were taken from 
the calcos pipeline “xld” hies, which include a summation over the four FP-POS exposures, 
eliminating bad pixels and suppressing hxed pattern noise such as grid-wire shadows. The 
visits 3-5 spectra for each grating setting were registered to visit 3 by cross-correlation, and 
co-added to increase S/N. In all cases, extracted lamp-hash lines from the individual FP- 
POS exposures were very stable in velocity, demonstrating the repeatability of the calcos 
pipeline wavelength assignments (although not necessarily their accuracy). 

Figure 8 compares selected prohles of bright lines from the various COS grating/detector- 
segment combinations against the STIS E140M spectrum from visit 1, convolved with the 
COS line-spread function. The zero point of the STIS wavelength scale (as in Fig. 2) was 
derived from the average velocity of O I A1304, A1306, and Cl I A1351, the brightest low- 
excitation narrow chromospheric lines recorded by STIS. The ground-state O I A1302 reso¬ 
nance line was avoided owing to a blue asymmetry imposed by redward interstellar atomic 
oxygen absorption (see Fig. 9). Note also that the STIS Al I A1670 resonance line is trun¬ 
cated in its red wing by a gap between adjacent echelle orders. The derived chromospheric 
radial velocity was —23 km s“^, compared with the photospheric value of —20 km s“^ ex¬ 
pected in 2012 near orbital phase 0.5 (Konig et al. 2005). This indicates that part of the 
velocity correction, amounting to —3 km s“^, could be a zero-point error in the STIS wave¬ 
length scale for this observation, assuming that the FUV chromospheric velocity measured 
by STIS, and the optical photospheric velocities measured by Konig et ah, refer to the same 
frame, and further that the optical measurements are freer of systematics. Aside from the 
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possible zero-point correction, the STIS wavelength scales have high internal precision. 

The COS spectra are depicted in two ways: (1) with the original pipeline wavelengths, 
adjusted only for the assumed stellar radial velocity (—20 km s“^); and (2) for a simple 
modihcation of the COS wavelength scales to force better global agreement with STIS. The 
latter corrections, expressed as equivalent velocities, are of the form: 

Av = VQ + a{X- Ao) , (1) 

The value Vq = —3 km s“^ was used for all the grating/segment combinations. The a 
coefficients and reference wavelengths, Aq, are as follows: —0.112 km s“^ and 1340 A 
for G130M/A; —0.054 km s“^ A“^ and 1400 A for G160M/B; and —0.046 km s“^ A“^ and 
1640 A for G160M/A. Here, the negative a coefficients are compensating for systematic 
redshifts of the assigned wavelengths that become progressively larger with wavelength. For 
example, the Si IV A1393 feature recorded by G130M/A (as in the earlier Ayres & France 
COS SNAPshot), would have a total instrumental redshift of 10 km s“^, about half of the 
Si IV velocities (~ 20 km s“^) measured in that study. However, a direct comparison between 
the two is subject to unknown systematic errors, not the least due to the different acquisition 
techniques, so the difference should not be taken too literally. 

The corrected COS spectra, in general, agree well with ground-truth STIS, although 
there is at least one important exception - Si IV A1393 at the extreme blue edge of G160M/B 
- where the linear velocity model appears to fall short. The cited grating-dependent velocity 
corrections were applied to the time-resolved COS prohles illustrated earlier in Figs. 4 and 5. 
Nevertheless, the fundamental velocity measurements of the brightest lines, at least for the 
quiescent periods of EK Dra, were derived from the more precise STIS spectra, as described 
below. 


3.2. Interstellar Features 

As mentioned earlier, there are interstellar absorptions present in several of the ground- 
state resonance lines of the STIS FUV spectrum of EK Dra. These offer the opportunity 
to check the STIS zero-point velocity deduced from narrow chromospheric emissions. The 
ISM features are redshifted with respect to EK Dra, with an average heliocentric velocity 
of about —2 km s“^. According to the maps of Redheld & Linsky (2008), EK Dra falls 
at the edge of the Local Interstellar Cloud (LIC) and also at the edge of the smaller North 
Galactic Pole (NGP) cloud. Depending on the velocity vector assumed for the LIC (e.g.. Cry 
& Jenkins 2014), the predicted heliocentric velocity in the direction of EK Dra falls in the 
range —(1-3) km s“^; while for the NGP cloud, the projected velocity is close to 0 km s“^: 
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Fig. 8.— Selected intervals from the various COS grating/detector-segment combinations, averaged 
over quiescent visits 3-5, compared to the STIS visit 1 spectrum, smoothed to COS resolution. The 
x-axis scale in each panel is AA (A) relative to the wavelength cited in the axis title. The STIS 
spectra, in the stellar frame defined by narrow chromospheric emission lines, are represented by 
red curves and open circles; the red dotted curves are the smoothed 1 a photometric noise (per 
resel). Blue dashed traces are for the uncorrected COS pipeline wavelengths, but accounting for 
the predicted photospheric radial velocity. Black dotted curves are for a simple modification of the 
COS wavelength scales to force better global agreement with STIS. Black dotted-dashed curves are 
smoothed COS 1 a photometric errors. For the STIS and COS spectra, only every fifth point is 
plotted. The narrow spike in the C I AI481 panel is a STIS hot pixel. 
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the galactic coordinates of EK Dra are nearly perpendicular to both directions of motion. It 
also is possible, and perhaps likely, that both clouds are contributing to the ISM absorptions 
in the EK Dra spectrum. 

Weighing in on the matter is a recent STIS NUV high-resolution echelle (E230H-2713) 
spectrum of EK Dra, taken in the 2014 ISM SNAP program mentioned earlier: a 600 s 
exposure through the 0.2"x0.2" slot, without a peak-up (see Table 1). The spectrum covers 
the important Mg II resonance lines at 2796 A and 2803 A, which have ISM absorptions on 
the red sides of their broad chromospheric emission cores. The Mg II ISM features fall at a 
velocity of —2 km s“^ heliocentric. However, cross-correlation against the same E230H-2713 
bandpass from a Cen A suggests that there might be a —1 km s“^ systematic error in the 
EK Dra NUV photospheric spectrum, which would revise the ISM velocity to about —1 km 
s“^. The FUV and NUV interstellar features are illustrated in Figure 9. The agreement 
between the two independent sets of ISM velocities provides additional conhdence in the 
zero-point velocity calibration of the STIS FUV spectra. 


3.3. STIS and COS Spectral Measurements 

Tables 3-5 list representative measurements from STIS visit 1 (Table 3), the COS visits 
3-5 averages for both G130M and G160M (Table 4), and the COS visit 2 flare (Table 5: 
G160M only). The STIS entries include important bright lines below the shortwavelength 
cutoff of COS G130M/A at 1290 A, such as H I Lya, Si III A1206, and the N V resonance 
doublet at 1240 A. STIS also recorded detailed prohles of the O I 1305 A triplet lines, which 
fall at the blue end of G130M/A but are corrupted by atomic oxygen skyglow in most of 
the COS pointings. The COS Tables include weaker features that were not detectable in 
the STIS spectrum, especially the key density-diagnostic multiplet of O IV] in the 1400 A 
region. The same G160M features in Table 4 also were measured in the integrated G160M 
flare spectrum of visit 2, and reported in Table 5, to compare with the more quiescent 
intervals following the outburst. 

For GOS, the measured Doppler shifts are affected by residual errors in the simple 
correction of the G130M and G160M wavelength scales, so should be treated with caution. 
However, the absolute fluxes and derived line widths, say of narrow components (NG) and 
broad components (BG) in a bimodal £t, should be more reliable, since neither depends 
sensitively on accurate wavelengths. Any systematic Doppler shifts between NG and BG 
within a given spectral feature should be reliable as well, because the (corrected) GOS 
wavelengths have sufficient precision over small wavelength separations. 
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Fig. 9.— Montage of ISM absorptions in the cores of strong chromospheric emissions of EK Dra, 
from the STIS FUV and NUV spectra. Each profile was normalized to the local flux density 
maximum in the velocity range, and separated vertically for clarity. The hatched emission feature 
in the center of the Lya absorption core is H I skyglow. The dotted-dashed line marks the predicted 
stellar radial velocity (not relevant for D I because H I dominates the emission profile), which was 
the same in the two separate epochs (FUV in 2012, NUV in 2014). Dotted lines delimit the range 
of predicted heliocentric velocities for the LISM and NGP clouds in the direction of EK Dra, which 
happens to fall at the edges of both structures, and nearly orthogonal to the respective velocity 
vectors. The NUV Mg II components are narrower than the FUV counterparts because the STIS 
high-resolution echelle was used at the longer wavelengths, and the thermal broadening is less for 
the heavier ion. 
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Table 3. HST Measurements: STIS FUV and NUV 


Species 

'^lab 

Flag 

Vh 

FWHMl 

/l 

fc 


(A) 


(km 

s-A 

(10“^'^ cgs) 

(lO-i-* cgs A-i) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

cm 

1175m 

INT 

(1175.7 A) 

(2.2 A) 

6.1 

0.00 

Si III 

1206.500 

GAU2 

+6±1 

69±5 

4.3±0.5 

0.00 

Si III 

1206.500 

GAU2 

+17±8 

184±33 

2.8±0.5 


D I 

1215.339 

GAUIS 

+ 18±0 

13±1 

-3.3±0.3 

67, -44 

HI 

1215.670 

GAU2’^ 

-8±1 

(231±7) 

(71±2) 

0.43 

HI 

1215.670 

GAU2’^ 

-4±4 

(565±21) 

(32±2) 


HI 

1215.670 

GAUV 

-9.3±0.7 

(272±3) 

(94±1) 

0.29 

HI 

1215.670 

GAUIX 

+20±0 




HI 

1215.670 

INT 

(1215.6 A) 

(5.0 A) 

67 

0.43 

N V 

1238.810 

GAU2D 

+4±1 

47±5 

0.81±0.12 

0.00 

N V 

1238.810 

GAU2D 

+8±6 

157±25 

0.97±0.11 

[0.477, -2.4] 

01 

1302.168 

GAU2R 

+3.9±0.7 

36±1 

1.6±0.1 

0.00 

01 

1302.168 

GAU2R 

+ 16±1 

15±1 

-0.35±0.12 


01 

1304.858 

GAUl 

-0.5±0.5 

33±1 

1.2±0.0 

0.00 

01 

1306.029 

GAUl 

+0.0±0.5 

31±1 

1.2±0.0 

0.00 

CII 

1334.532 

GAU2R 

+6.1±0.4 

51±1 

4.5±0.1 

0.02 

CII 

1334.532 

GAU2R 

+ 17±0 

14±1 

-1.2±0.2 


CII 

1335.702 

GAU2 

+3.0±0.4 

48±2 

3.7±0.2 

0.02 

CII 

1335.702 

GAU2 

+ 11±2 

137±8 

2.9±0.2 


CII 

1334m 

INT 

(1335.2 A) 

(2.8 A) 

11 

0.02 

CIU 

1351.656 

GAUl 

+0.4±0.9 

20±2 

0.23±0.02 

0.03 

FeXXI 

1354.08 

GAUl 

-4±3 

103±8 

0.68±0.04 

0.03 

01 

1355.598 

GAUl 

-1±2 

27±4 

0.13±0.01 

0.03 

Cl 

1355.844 

GAUl 

0±2 

27±4 

O.lOiO.Ol 

0.00 

Si IV 

1393.760 

GAU2D 

+6.2±0.6 

51±2 

2.7±0.1 

0.02 

Si IV 

1393.760 

GAU2D 

+ 10±3 

174±8 

2.6±0.1 

[0.546,-0.9] 

CIV 

1548.204 

GAU2D 

+7.1±0.6 

50±2 

5.6±0.3 

0.03 

CIV 

1548.204 

GAU2D 

+ 11±2 

179±8 

6.9±0.3 

[0.585,-2.3] 

He II 

1640.404 

GAU2 

+5.8±0.7 

54±3 

8.3±0.7 

0.11 

He II 

1640.404 

GAU2 

+28±14 

132±36 

2.4±0.7 


STIS Visit 1: N V + Si IV + C IV 

Narrow component 

GAU2 

+5.1±0.1 

50±1 

0.48±0.01 
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Table 3—Continued 


Species 

-^lab 

Flag 

Vh 

FWHMl 

/l 

fc 


(A) 


(km s 

-b 

(10-14 cgs) 

(10-14 cgs A-i) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Broad component 

GAU2 

+8.7±0.4 

137±1 

0.52±0.01 



Mg II Chromospheric Emission Cores in 

Stellar Rest Frame 

Mg II 

2796.352 

GAUl 

+0.1±0.3 

60.0±0.9 

150±2 

13.2 

Mg II 

2803.531 

GAUl 

+0.3±0.3 

54.3±0.6 

117±2 

13.2 



Mg II Interstellar Lines 

in Heliocentric Frame 


RSI 

2796.352 

GAUIS 

-1.8±0.1 

5.9±0.3 

-13.8±0.6 

201,-531 

III 

2803.531 

GAUIS 

-1.2±0.1 

5.4±0.3 

-8.6±0.4 

166, -499 


Note. — Cols. 1 and 2 list the species and laboratory (vacuum) wavelength: suffix “m” 
indicates a multiplet. Col. 3 is a processing flag: see text for specific descriptions. Cols. 4 
and 5 are centroid and width of the Gaussian profile, in velocity units. Col. 6 is the 
integrated flux, in units of 10“^^ erg cm“^ s“^. Col. 7 usually is the derived continuum 
flux density, with units of 10“^^ erg cm“^ s“^ A“^. If two values are listed, the second 
value is the gradient of a sloping fit (see text). Note, in the specific case of a “GAU2D” 
constrained doublet fit, the first set of entries are the continuum parameters, but the 
second, bracketed, set is the ratio of the integrated flux of the secondary member divided 
by the principal member, and the deviation of the measured wavelength difference relative 
to the laboratory value (in km s“^). An error of zero indicates an uncertainty less than 
the least significant figure of the measurement. 

^'Fit ignores deep ISM central reversal, solely to obtain centroid velocity of outer profile. 












Table 4. HST Measurements: COS Visits 3-5 (Quiescent) 


Species 

-^lab 

Flag 

Vh 

FWHMl 

/l 

/c 


(A) 


(km s 

-A 

(10-14 cgs) 

(10-14 cgs A-i) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 


G130M 

CII 

1334.532 

GAU2R 

-k9.7±0.0 

60±0 

4.5±0.0 

0.03 

CII 

1334.532 

GAU2R 

-k22±0 

17±0 

-0.97±0.04 


CII 

1335.702 

GAU2 

-k0.2±0.2 

48±0 

3.2±0.0 

0.03 

CII 

1335.702 

GAU2 

-kl6±0 

122±1 

3.4±0.1 


CII 

1334m 

INT 

(1335.2 A) 

(3.2 A) 

11 

0.03 

Cll 

1351.656 

GAUl 

-k0.7±0.2 

26±1 

0.23±0.00 

0.04 

FeXXI 

1354.08 

GAUl 

-k5.6±0.7 

120±2 

0.52±0.01 

0.04 

01 

1355.598 

GAUl 

-2.1±0.7 

38±1 

0.15±0.00 

0.04 

Cl 

1355.844 

GAUl 

-k2.1±0.7 

38±1 

O.lliO.OO 

0.00 

0 V 

1371.296 

GAUl 

-k6.0±0.9 

53±2 

0.12±0.00 

0.05 

Si IV 

1393.760 

GAU2D 

-k4.7±0.2 

47±0 

2.1±0.0 

0.06 

Si IV 

1393.760 

GAU2D 

-kl4±0 

133±1 

2.9±0.0 

[0.5248,-1.1] 

OIV 

1401.157 

GAU7 

+9±1 

91±2 

0.20±0.00 

0.07 

OIV 

1397.202 

GAU7 



0.008±0.003 


OIV 

1399.769 

GAU7 



0.065±0.004 


OIV 

1404.780 

GAU7 



0.052±0.004 


SIV 

1404.808 

GAU7 



0.014±0.001 


SIV 

1406.051 

GAU7 



0.046±0.004 


OIV 

1407.379 

GAU7 



0.064±0.004 


G160M 

Si IV 

1393.760 

GAU2D 

-1.7±0.2 

60±1 

2.5±0.0 

0.06 

Si IV 

1393.760 

GAU2D 

-k4.9±0.4 

167±2 

2.6±0.0 

[0.552,-k7.1] 

OIV 

1401.157 

GAU7 

-klOil 

92±3 

0.21±0.01 

0.07 

OIV 

1397.202 

GAU7 



0.008±0.003 


OIV 

1399.769 

GAU7 



0.055±0.004 


OIV 

1404.780 

GAU7 



0.038±0.004 


SIV 

1404.808 

GAU7 



0.014±0.001 


SIV 

1406.051 

GAU7 



0.046±0.004 


OIV 

1407.379 

GAU7 



0.063±0.005 


Si II 

1526.707 

GAUl 

-k0.2±0.4 

50±1 

0.60±0.01 

0.10 



























Table 4—Continued 


Species 

'^lab 

Flag 

V-L 

FWHMl 

/l 

/c 


(A) 


(km s' 


(10“^"^ cgs) 

(10-14 cgs A-i) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Si II 

1533.431 

GAUl 

+3.2±0.4 

51±1 

0.67±0.01 

0.10 

CIV 

1548.204 

GAU2D 

+2.6±0.2 

48±0 

4.2±0.1 

0.14 

CIV 

1548.204 

GAU2D 

+ 18±0 

151±1 

6.8±0.1 

[0.564,+3.3] 

He II 

1640.404 

GAU2 

+4.1±0.2 

50±1 

5.9±0.1 

0.15 

He II 

1640.404 

GAU2 

+20±1 

128±2 

4.5±0.1 


Cl 

1657m 

INT 

(1657.1 A) 

(3.0 A) 

5.8 

0.21 

Alll 

1670.787 

GAU2 

+ 1.9±0.7 

33±3 

0.52±0.07 

0.18 

Alll 

1670.787 

GAU2 

-1±1 

107±7 

l.OiO.l 



Note. — See comments table 3. 







Table 5. HST Measurements: COS Visit 2 (Flare) 


Species 

Alab 

Flag 

VL 

FWHMl 

/l 

/c 


(A) 


(km s 


(lO-i-* cgs) 

(lO-i-* cgs A-i) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Si IV 

1393.760 

GAU2D 

+0.8±0.6 

61i2 

9.1i0.3 

0.39 

Si IV 

1393.760 

GAU2D 

+31±2 

207i4 

9.6i0.3 

[0.562, i6.2] 

OIV 

1401.157 

GAU7 

+19±3 

124i8 

l.liO.l 

0.39 

OIV 

1397.202 

GAU7 



0.07i0.03 


OIV 

1399.769 

GAU7 



0.45i0.04 


OIV 

1404.780 

GAU7 



0.25i0.04 


SIV 

1404.808 

GAU7 



0.06i0.01 


SIV 

1406.051 

GAU7 



0.20i0.04 


OIV 

1407.379 

GAU7 



0.42i0.04 


Si II 

1526.707 

GAUl 

Oil 

53i3 

1.5i0.1 

0.53 

Si II 

1533.431 

GAUl 

ili2 

58i4 

1.6i0.1 

0.53 

CIV 

1548.204 

GAU2D 

i5.7i0.4 

45il 

26i0 

0.58 

CIV 

1548.204 

GAU2D 

i41il 

197i2 

43i0 

[0.583,-0.2] 

He II 

1640.404 

GAU2 

i6.8i0.4 

50il 

26il 

0.69 

He II 

1640.404 

GAU2 

i36i2 

164i5 

18il 


Cl 

1657m 

INT 

(1657.1 A) 

(3.0 A) 

13 

0.74 

A/II 

1670.787 

GAU2 

-2il 

27i4 

1.6i0.2 

0.77 

Alll 

1670.787 

GAU2 

il9i5 

169il4 

4.4i0.3 



Note. — See comments Table 3. 
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Reference (“laboratory”) wavelengths (in vacnnm) were taken from Ayres (2015), except 
for the solar flare line [Fe XXI], which was not inclnded in the a Cen stndy. The NIST 
Atomic Spectra Database (ASD; Kramida et ah 2013)0 wavelength (1354.08 A) was adopted 
for this key coronal forbidden transition, bnt the qnality of the measnrement is rated poor. 
Recently, Yonng et ah (2015) proposed a revised wavelength of 1354.106±0.023 A, based on 
measurements of [Fe XXI] in solar flares. If correct, it would shift the [Fe XXI] velocities in 
Tables 3-5 by 6±5 km s“^ to the blue. 

A variety of fitting strategies were employed, following the lead of the Ayres (2015) 
study of STIS FUV echelle spectra of a Cen A and B. The particular measurement scenario 
is flagged explicitly in col. 3 of the Tables, as follows. “INT” refers to a numerical flux 
integration over the velocity bandpass listed in col. 5. The value in col. 7 is a derived 
continuum flux density, assumed constant over the spectral feature or multiplet, based on 
filtering the surrounding intensities to suppress emission lines. “GAUl” refers to a single 
Gaussian fit. “GAUIS” is a special case in which a sloping continuum was allowed: col. 7 
now contains the zero point and gradient of a linear continuum fit: (/c)a = Gq -|-Gi (A — Al), 
where fc is the monochromatic continuum flux density and Al is the transition wavelength. 
“GAUIX” is another specialized fit, applied to the ISM absorption core of H I Lya. The 
fitting function is of the form where x = AA/AAd and AAd is the “Doppler width.” 
Gompared with the normal Gaussian (e“*^), the steep-sided flat-bottomed GAUIX prohle 
is more appropriate for the saturated Lya absorption core. However, the FWHM and total 
flux have little meaning, especially because the imposed continuum level (at the apparent 
peak of Lya) is rather arbitrary, so the derived values are not reported in the table; only the 
velocity shift is significant. 

“GAU2” is a bimodal Gaussian fit of the narrow/broad variety, used specifically for 
the N V -|- Si IV -|- G IV collective filtered profiles, but also for other isolated features such 
as Si III A1206. For the filtered hot-line average, the component fluxes are narrow/total 
and broad/total. “GAU2R” refers to a double-Gaussian fit applied to ostensibly unrelated 
features (say, stellar emission and superimposed interstellar absorption), for which the second 
component is negative in flux to mimic a central reversal or ISM absorption. In these cases, 
the width of the principal component should not be taken too literally, although the velocity 
shift should be reliable. 

“GAU7” is a 7-component Gaussian model covering combined multiplets of similar 
transitions, namely O IV] -|- S IV], for which the widths of the features were constrained to 
be the same, and the Doppler shifts as well (i.e., the wavelength differences between the 


^see: http://www.physics.nist.gov/PhysReilI)ata/ASD/lines_form.htinl 
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components were fixed at the reference values, but a shift of the whole transition array was 
allowed). The derived collective velocity shift and average prohle width, together with the 
associated uncertainties (and the continuum flux) are listed for the principal transition. For 
the remaining multiplet members, only the unique flux values are provided. 

Finally, “GAU2D” refers to a specialized constrained bimodal Gaussian £t applied to 
the hot-line doublets, such as Si IV A1393 -|- A1402. Figure 10 illustrates such fits to the (un¬ 
smoothed) EK Dra hot lines from the co-added STIS spectrum. This bimodal strategy for 
decomposing the complex hot-line profile shapes of cool stars has a long history, as discussed 
recently by Ayres (2015) in the context of low-activity a Gen A (and companion B [K1 V]), 
and for the earlier GOS pointing on EK Dra itself by Ayres & France (2010). The upper 
paired panels of Fig. 10 illustrate the application to the N V, Si IV, and G IV resonance 
doublets: principal components to the left, secondary to the right (but with one-half the 
y-axis scale). Under optically-thin conditions, the weaker components would be half the flux 
densities of the stronger ones, and thus should appear similar in these scaled diagrams. That 
expectation is met for N V, but becomes progressively out of sync for Si IV and G IV. The 
fitting procedure assumed that the doublet members were self-similar: NG of the same width, 
BG of the same width; Doppler shifts of the NG the same, and similarly for the BG (although 
possibly different between narrow and broad); and the same flux ratios, narrow/broad. Only 
the relative total fluxes between the doublet members was allowed to float, as well as the 
wavelength difference between the two transitions, to allow for small errors in the laboratory 
values, or subtle residual distortions in the STIS wavelength scales. 

For this GAU2D case, the first set of parameters listed (in Tables 3-5) is for the NG 
of the constrained doublet fit, and col. 7 holds the continuum flux as usual. The col. 2 
wavelength is for the principal member of the doublet. The second set of parameters listed 
is for the BG. However, highlighted by square brackets, col. 7 now has the doublet flux ratio, 
weaker/stronger, which should be close to 0.5 in the optically thin limit (for the specific 
hot-line doublets here); as well as the deviation of the fitted wavelength difference from the 
reference value, expressed in km s“^. These latter shifts are small for the STIS spectrum, 
but can be substantial in several of the GOS cases, mainly reflecting residual errors in the 
simple correction of the GOS wavelengths. 

The bottom panel of Fig. 10 illustrates a GAU2 bimodal fit to a filtered sum of the 
three hot-line doublets, each member (six altogether) scaled to the same integrated flux over 
the interval |u| <65 km s“^. The scaled fluxes were combined on a common velocity grid, 
applying an Olympic filter (throwing out the highest and lowest of the six values) for the 
central ±65 km s“^ zone, and a more aggressive filtering (eliminating the three highest values) 
for the outer intervals where the influence of weak emission blends on the broad component 
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VELOCITY (km s"’) 


Fig. 10.— Bimodal Gaussian fits to the EK Dra hot-line doublets, and a filtered average of 
the three, from the co-added STIS spectrum. In the upper paired panels, black dotted-dashed 
curves are smoothed 1 a photometric errors (per resel). In the bottom panel, the black dotted 
dashed curves indicate ±3 standard errors of the mean (s.e.) of the filtered average (per bin). The 
filtering was designed to suppress accidental emission blends from unrelated species. The small 
s.e.’s demonstrate that the profiles are close to self-similar, aside from the occasional blend. Open 
circles represent the STIS spectra: only every other point is plotted. In the upper panels, the two 
members of each doublet were fitted simultaneously, imposing various constraints on the narrow 
and broad sub-component widths, shifts, and strengths (see text). The NC are blue, the BC are 
red, and the sum is green. The STIS E140M line spread function for the photometric aperture 
(FWHM~ 8 km s“^: yellow shaded prohle in the bottom panel) was accounted in the modeling 
procedure. Note the factor of 2 y-axis scale change between the upper left panels, for the principal 
doublet members, and the right panels, for the secondary members. With that scaling, the two 
features should have similar appearance in the optically thin limit. 
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is strongest. As with the individual bimodal hts, the hybrid prohle has a relatively narrow 
NC, about 50 km s“^ FWHM, and is redshifted by 5 km s“^, similar to those of low-activity 
G dwarfs like a Cen A; but a very broad BC, about three times wider than the NC, and 
more redshifted (9 km s“^). The enhanced broadening and redshift of EK Dra’s BC are the 
main distinguishing characteristics compared with low-activity counterparts like a Cen A 
(e.g., Fig. 2). 

Uncertainties on the htted parameters were estimated by a Monte Carlo approach: a 
reference prohle (the derived ht itself) was perturbed by random realizations of the smoothed 
photometric error, and then rehtted. The standard deviations of the derived parameters over 
many such trials were taken as a fair gauge of the measurement uncertainties devolving from 
the photometric noise component. 


3.4. The Isolated FUV Continuum Bursts 

Recall the two bursts (at t =3.6 hr and 3.7 hr) in the G160M blue and red continuum 
bands, which remarkably were not mirrored in any of the other activity tracers. These 
isolated bursts are distinctly bluer (indicating a hotter continuum) than the more equal count 
rates of the two separated continuum bands later on, in quiet visits 3-5. This continuum- 
only behavior is challenging to understand, and seems to have few, if any, precedents among 
existing records of stellar flare events. This likely is because there are only a handful of stellar 
flare campaigns simultaneously at UV and soft X-ray energies. One example is Mitra-Kraev 
et al. (2005), who were able to use the Optical Monitor on XMM-Newton, with broad-band 
FUV and NUV hlters, to jointly record a number of UV/X-ray events on several dMe stars. 
The authors found that the UV rise typically preceded the soft X-rays by several minutes, 
which usually is the case for white-light flares (WLF) on the Sun. There is perhaps only a 
single example in their several time series where a UV burst was seen during a flare decay 
without an X-ray counterpart (their Fig. 1). On the Sun, WLFs are difficult to separate from 
the glare of the photosphere, but if one is seen, it usually will have an extensive supporting 
context, especially at high energies. Fang & Ding (1995) divided solar WLFs into two classes: 
Type I, or “prompt,” events which are coincident with the hard X-ray (HXR) burst, and 
apparently result from direct excitation of the chromospheric gas by the flare’s energetic 
particle beams; and the less frequent Type II, or “delayed,” bursts, which occur several 
minutes after the HXR impulsive phase, and appear to be deeper seated, in the photosphere 
itself. Fang & Ding did not speculate on a source for the Type IIs. Perhaps the delayed 
FUV continuum brightenings on EK Dra are related to the solar WLF Type IIs. One idea, 
motivated by the hot-line redshifts during the flare decay, is that the bursts are a response to 
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the ballistic splash-down of the cooling gas onto the photospheric layers; something analogous 
to accretion streams in T-Tauri stars. 


3.5. The [Fe XXI] Coronal Forbidden Line 


Ayres et ah (2003) carried out a search for coronal forbidden lines in STIS FUV spectra 
of the late-type objects with sufficiently deep exposures available at the time. In the hyper¬ 
active members of the sample, like RS Canum Venaticorum binaries and dMe flare stars, 
the intrinsically faint [Fe XXI] A1354 feature often was visible. Furthermore, the [Fe XXI] 
line strength was found to be linearly correlated with the ROSAT soft X-ray flux. In the 
moderate-to-slow rotators (in terms of nsini), the average FWHM of the coronal iron line 
was about 110 km s“^, somewhat larger than the expected thermal width of 90 km s ^ at 
the formation temperature of 10 MK. In several of the fast rotators, there was a suggestion of 
~ 2x enhanced broadening beyond that expected from nsini alone. The “super-rotational” 
broadening was attributed to highly extended hot coronal structures. 


In EK Dra, the COS visits 3-5 average prohle of [Fe XXI] does not display obviously 
extended red wings, like the lower temperature subcoronal features, although certainly the 
coronal forbidden line is very broad, approximately mid-way in width between the hot-line 
NC and BC. Nevertheless, the FWHM~ 120 km s“^ is close enough to the low-n sini average 
noted above that it would be difficult to make a case for super-rotational broadening. 

The COS velocity of [Fe XXI] in the visits 3-5 average spectrum is -|-6 km s“^ for the 
NIST reference wavelength, although with the Young et al. (2015) revised value, the shift 
would decrease to about 0 (±5 km s“^ from the ax of the reference wavelength). At the same 
time, the STIS velocity of [Fe XXI] would shift to —10 km s“^ (±6 km s“^, where the larger 
a partly is due to the reference wavelength, and partly to the low S/N of the measurement). 

Because [Fe XXI] correlates well with the ROSAT flux for hyper-active coronal sources, 
the measured intensity during the pre- and post-flare intervals can be exploited to estimate 
the quiescent X-ray flux of EK Dra during the 2012 campaign, noting that the existing 
ROSAT survey Lx = 10^®'® erg s“^ (0.1-2.4 keV) reported by Giidel et al. (1995) is not 
exactly simultaneous. The scaling relation of Ayres et al. (2003) can be expressed as: 


2x10'' 


. /peXXI. 
/bol 


^bol 


( 2 ) 


The stellar bolometric flux, /boi ~ 2.45x10“^ erg cm“^ s“^, can be estimated from V and a 
bolometric correction derived from B — V (see Ayres et al. 2005). The bolometric luminosity, 
Tboi ~ 3.4x 10^^ erg s“^, follows from the distance (34 pc), and is similar to Lq. The [Fe XXI] 


rsj 
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flux ratio, based on the STIS/COS average /l ~ 6x10“^^ erg cm“^ s“^, is 2.5x10“’^: in the 
upper tier among the coronal objects considered by Ayres et ah (2003). Finally, the inferred 
X-ray luminosity is Lx ~ 1.7x10^° erg s“^, about twice the ROSAT survey value, but 
certainly consistent within the uncertainties of the [Fe XXI] scaling law, and the temporal 
variability of Lx noted by Giidel et ah (1995). 


3.6. Correlations among Bright Lines in the Flare and Quiescent Periods 

Figure 11 is a correlation diagram pitting various diagnostics against the Si IV dou¬ 
blet, which was recorded in every COS observation. As with the previous diagram(s), the 
individual flux values from all the visits were normalized to the corresponding averages over 
visits 3-5. The upper distributions of symbols are contributed by the large flare, while the 
clump at (1,1) is dominated by the quiescent intervals, which nevertheless display low-level 
variability including a small flare decay. The straight lines are power-law hts to the indi¬ 
vidual distributions. The derived slopes range from 1.44 and 1.38 for [Fe XXI] and C IV, 
respectively, to slightly above unity for He II (1.10), and slightly below for chromospheric 
C II (0.89). These slopes mirror similar power laws seen between, say. X-rays and C IV 
in diverse samples of solar-type dwarfs (e.g., Ayres et ah 1995), and in this instance must 
reflect the time-evolution of the post-flare cooling process. 


3.7. Quiescent and Flare Profiles of C IV 

Figure 12 compares C IV doublet prohles from the COS C160M exposures of the visit 2 
flare decay, the visits 3-5 average COS G160M spectrum, and a multi-epoch STIS average 
for the quiet solar twin a Cen A (Ayres 2015); all including bimodal constrained Gaussian 
hts. The quiescent prohles of EK Dra show redshifts of both components, less for the NC 
(~ -|-6 km s“^, for the average of the STIS C IV doublet members), more for the BC 
(~ -|-10 km s“^, for the STIS doublet average). This compares with a somewhat smaller 
redshift of both a Cen A components (~ -|-5 km s“^ for each). The width of the EK Dra 
C IV NC in the quiet spectrum (50 km s“^) is somewhat larger than in a Cen A (36 km 
s“^), but the BC is much broader (180 km s“^ vs. 66 km s“^). In the hare, remarkably, the 
C IV NC has a similar redshift and FWHM (-1-6 km s“^ and 45 km s“^, respectively) to the 
“quiet” average, and indeed to the low-activity solar twin a Cen A. At the same time, the 
hare BC not only is broader than in the quiet cases (~ 200 km s“^), but is substantially 
more redshifted (at about -|-40 km s“^). The larger BC redshifts during the hare are notable 
but not unusual; similar behavior, with similar velocity shifts, has been seen previously in 
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Fig. 11.— Correlation diagram pitting various diagnostics from all the visits, normalized to the 
visits 3-5 average, against the similarly normalized Si IV doublet count rates. The systematic trends 
toward the upper right represent the flare decay in visit 2, while the dispersions of the clusters of 
points near the origin (1,1) indicate the typical variability during quieter times. 
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STIS spectra of flares on the dMe stars AU Mic (Robinson et al. 2001) and AD Leo (Hawley 
et al. 2003), and appears to be a signature of the post-flare cooling process. What perhaps 
is more remarkable is that despite the very large difference in emission levels of the two solar 
analogs (expressed, say, in terms of surface fluxes), the “quiet” C IV prohles are as similar 
as they are, not only in the bimodal nature of the line shape asymmetries, but also in the 
magnitudes of the velocity shifts and line widths. 


3.8. Densities from the O IV] Semi-Forbidden Mnltiplet 

Flower & Nussbaumer (1975) originally proposed that flux ratios within the O IV] 
1400 A intercombination mnltiplet are sensitive to density over the range rie ~ 10^° — 10^^ 
cm“^ thought to be characteristic of solar conditions - quiescent at the low end, flaring on 
the high side - at the ionization equilibrium peak temperature of ~ 10^ K. There is 
a long history of applications to solar and stellar plasmas (see, e.g.. Cook et al. 1995, and 
references therein), continuing to the present (e.g., Keenan et al. 2009; Ayres 2015). Figure 13 
illustrates the O IV] region from the initial visit 2 flare interval (G160M, only); the quiet 
averages (visits 3-5) from both G160M and G130M; and the a Cen A STIS spectrum. The 
logarithmic presentation emphasizes the wider, more redshifted broad components of the 
EK Dra prohles compared with low-activity a Cen. The 7-component Gaussian hts (hve 
O IV], two S IV] transitions) take into account the small contamination of O IV] A1404 by 
S IV] A1404, scaled from S IV] A1406 (one of the htted features) according to plasma emission 
models (see Ayres 2015). 

Table 6 summarizes O IV] diagnostic ratios constructed from the COS measurements 
in Tables 4 and 5, and including values from the multi-epoch average a Cen A STIS spec¬ 
trum. The nomenclature for the ratios is that of Ayres (2015): the numerators in the 
traditional density-sensitive i?i = A1407/A1401 and i ?2 = A1407/A1404 were replaced by 
|(A1399 -|- A1407) to improve S/N with these generally faint emissions: A1399 and A1407 
share a common upper level and their ratio is very close to unity according to theory. The 
hybrid numerator ratios are designated R*. The ratio A1397/A1404 also is composed of 
transitions that share a common upper level, with a theoretical value of 0.13 (Keenan et al. 
2009). However, in this case there is no advantage to form a hybrid flux between A1397 and 
A1404 because the former transition usually is too low in S/N, if detected at all, in normal 
astrophysical sources. 

The a Cen A1399/A1407 and A1397/A1404 ratios are close to the theoretical expecta¬ 
tions. This also is the case for EK Dra, although the uncertainties are larger, especially 
for marginally detected A1397. Note, however, that the diagnostic ratios Ri* and R 2 * are 
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Fig. 12.— C IV doublet profiles, and bimodal fits, from the EK Dra flare (COS visit 2: upper 
panel) and the quiet average (COS visits 3-5: middle panel), compared to the low-activity sunlike 
star a Cen A (from STIS: lower panel). Only every third point (black dots) is plotted for these 
oversampled broad profiles. Note the large scale change between the upper (flare) and middle 
(quiescent) panels for EK Dra. Black dotted-dashed curves are smoothed 10 a photometric errors 
(per resel). The EK Dra flare profiles are remarkably similar to the quiet average, despite the 
large flux increase, except for the more redshifted broad components. The C IV NC of EK Dra are 
somewhat similar to those of a Cen A, but the BC are much broader and occupy a larger fraction 
of the total doublet flux. 
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Fig. 13.— Density diagnostic O IV] (plus S IV]) 1400 A multiplet, together with multi-component 
Gaussian model (thick curves). Same arrangement of spectra as in Fig. 12. Note the logarithmic y- 
axes, and large scale change between the two EK Dra panels. As with the previous figure, only every 
third point (dots) is plotted. Dotted-dashed curves again are smoothed 10 d photometric errors. 
Note also that the middle panel contains both G130M (red squares for data; green curve for fit; red 
dotted-dashed curve for 10 d level) and G160M (black dots; blue dashed; black dotted-dashed) for 
the visits 3-5 average spectra, while the upper panel has only G160M (from the initial flare interval 
in visit 2). In the G130M/G160M comparison, the agreement between the wavelength-corrected 
prohles of Si IV A1402 is excellent, while the shorter wavelength component, A1393, displays a 
noticeable blue deviation in G160M, as noted earlier in Fig. 8. The numerous narrow lines of 
chromospheric species (e.g., S I, Fe II, and Ni II) seen in low-activity a Cen A are expected to be 
much less prominent in hyper-active EK Dra, owing to the general decline in importance of the 
lower temperature layers, relative to the hotter subcoronal plasma, with increasing activity. 











Table 6. O IV] Line Ratios 


Dataset 


i?2* 




4(A1399 +A1407) 

1(A1399 +A1407) 

A1399 

A1397 


A1401 

A1404 

A1407 

A1404 

(1) 

(2) 

(3) 

(4) 

(5) 

Theoretical Limits 

0.16-0.42 

0.3-2.5 

1.02 

0.13 

Electron Density Range (cm“^) 

lOlO _ 1012 

IQio _ 1012 

7^ 

Tie 

EK Draconis 

Quiet Average: COS G130M 

0.32±0.02 

1.2±0.1 

l.OiO.l 

0.15±0.05 

Quiet Average: COS G160M 

0.28±0.02 

1.6±0.2 

0.9±0.1 

0.21±0.08 

Flare: COS G160M 

0.40±0.04 

1.7±0.3 

l.liO.l 

0.27±0.11 

a Centauri A 

Quiet Average: STIS E140M 

0.185±0.002 

0.429±0.006 

0.99±0.03 

0.127±0.007 


Note. — Theoretical values from Keenan et al. (2009). 
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significantly elevated in EK Dra compared with a Cen. This is a clear indication of en¬ 
hanced electron densities in the more active star. The a Cen ratios fall near the low-density 
limits of the respective i?’s, indicating Ug just below 10^° cm“^; while the EK Dra ratios 
both are consistent with densities in the middle range, near or above 10^^ cm“^, an order 
of magnitude larger. Curiously, the inferred densities are not much different between the 
more quiet periods of the EK Dra COS campaign and the large FUV flare decay at the 
beginning; especially for R 2 *, which has the largest grasp of the two diagnostic ratios (i.e., 
largest variation in R over the range of density sensitivity). 

The O IV] densities inferred here are about a factor of 2 higher than those estimated 
by Gudel et ah (1995) from emission-measure modeling of the ROSAT survey pulse-height 
spectra, at least for the softer (2 MK) component of their 2-temperature £t (this is the com¬ 
ponent that showed rotational modulations, as mentioned earlier; whereas the harder, more 
dominant, ~ 10 MK one showed essentially no phase-correlated variability). Nevertheless, 
the agreement must be tempered by the fact that the subcoronal gas is only one-tenth the 
temperature of even the cooler coronal plasma component, and likely is not co-located. 

The high densities during the flare decay in visit 2 suggest short cooling times {t ~ 
(rze)”^), at least in the hot-line temperature range (~ 10^ K) where the radiative cooling 
curve peaks. This seemingly conflicts with the hours-long progression of the FUV decay. 
A similar dichotomy confronts interpretations of long-duration flares on the Sun and other 
late-type stars. The usual explanation is that the long duration events are analogous to 
two-ribbon (2-R) flares on the Sun (e.g., Giidel 2004), which involve a sequential cascade 
of excitation through large arcade-like magnetic structures, rather than an impulsive burst 
in a single coronal loop. Such 2-R events tend to be quite hot, exceeding 10 MK in peak 
temperature, with coronal electron densities up to about 10^^ cm“^ (e.g., Giidel 2004). 
Importantly, the decay of the flare plasma in a 2-R event can be dominated by conductive 
energy exchange with the lower cooler layers, rather than direct radiative losses at the super- 
hot coronal temperatures (see Aschwanden et al. 2008). Thus, the bright FUV emissions and 
subsonic plasma draining seen in the hot lines of EK Dra during the large FUV flare decay 
likely are a byproduct of these cooling processes playing out in the much more energetic 
setting of the hyper-active young solar analog. This brings to mind, as intimated by Ayres 
& France, the dynamic “coronal rain” phenomenon associated with post-flare loops on the 
Sun, as well as solar active regions in general (see, e.g., Antolin & Rouppe van der Voort 
2012, and references therein). 
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3.9. The EK Draconis Event in a Broader Context 

Because of the unique nature of the time-resolved FUV spectroscopy of the serendipitous 
EK Dra flare decay, it is important to place the event within the context of the nearly 
100 large stellar flares captured in the soft X-ray band by historical and contemporary 
high-energy observatories (as documented in the review by Giidel 2004). Aschwanden et 
ah (2008) have htted various combinations of the derived emission measures (EM), X-ray 
luminosities (Lx) at the flare temperature peak (Tp), and decay times (tf) with scaling 
relations approximated as power laws. Inverting their eq. 12 for the dependence of decay 
time on flare peak temperature yields an estimate of the (here unknown) Tp given the (known) 
tf. 

Tp ~ 10 MK(^^)°-® . (3) 

Given the apparent, albeit FUV, flare decay time of about 7 ks, the peak X-ray temperature 
would be about 60 MK. Using this value in eq. 13 of Aschwanden et ah yields Lx ~ 10^^'^ 
erg s“^, which is an order of magnitude larger than the “quiescent” Lx inferred from scaling 
the contemporaneous [Fe XXI] flux from the non-flare intervals. Supporting that estimate 
is the fact that the peak enhancement of G IV seen in GOS visit 2 is about a factor of 10, 
and the [Fe XXI] proxy for the high-temperature soft X-ray flux - although recorded later in 
visit 2, not simultaneously with G IV - parallels the G IV trajectory during the flare decay, 
suggesting that [Fe XXI] would have been enhanced by a similar large factor, and so too the 
soft X-ray flux. 

The inferred values of flare peak temperature. X-ray luminosity, decay time, and density 
(again with the caveat that the latter is from the cooler FUV plasma) place the EK Dra 
event in the upper-middle range of the documented stellar X-ray flares, as summarized by 
Aschwanden et ah; perhaps typical of a young active sunlike star, but more accessible in 
the case of EK Dra by its proximity: several times closer to Earth than the members of the 
similar-age young clusters Pleiades and Alpha Persei. 


4. CONCLUSIONS 

The joint STIS/GOS study of EK Dra was adversely impacted by the unfortuitously- 
timed large FUV flare (interrupting the intended instrumental cross-calibration), but never¬ 
theless the event itself provided valuable insight into the post-flare cooling process, mainly 
thanks to the high spectral resolution achieved for a variety of key chromospheric and higher 
temperature features. The apparent quiescent periods following the flare also contributed 
signihcant insight, not only with regard to the velocity cross-calibration of the two UV spec- 
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trographs, but also in the lack of phase-correlated prohle variations caused by rotational 
modulations. 

To summarize: (1) the FUV redshifts appear to be a characteristic of the normal surface 
activity of the star, rather than an accident of Doppler Imaging; (2) the hot-line prohles did 
not change qualitatively during the flare outburst (ten-times enhancement of C IV), except 
for the larger redshifts of the broad components; and (3) O IV]-derived densities were not 
very different between the quiescent periods and the FUV flare, despite the fact that the 
high densities themselves are characteristic of large stellar soft X-ray events (see, e.g., Benz 
& Giidel 2010). 

All of these properties suggest that the normal surface conditions on the star are not 
far removed from those in a big flare; in other words, the stellar corona might be dominated 
by flaring structures, with a distribution weighted toward the smaller events, as on the Sun, 
giving the appearance of quiescence at times, when exactly the opposite is true (see Giidel 
2004 for a review of the “flares all the way down” coronal scenario; what might be described 
as a “flare-ona”). 

One notable outcome of the time-resolved flare spectroscopy was the lack of strongly 
blueshifted emission, which naively might be anticipated from the catastrophic energy re¬ 
lease, and indeed such outbursts on the Sun are violently eruptive. However, as has been 
noted often in the flare literature, during the explosive heating of the lower atmosphere in 
an event, most of the outwardly expanding material undoubtedly is too hot to be visible 
in the FUV lines. Only the subsequent decay phase would contribute strongly to the FUV 
features, as the super-hot coronal plasma in the flaring arcades conductively redistributes 
heat into the lower atmosphere, where it is radiated away partially at FUV temperatures 
and inspires the subsonic plasma draining phenomenon manifested in the hot-line redshifts. 
Unfortunately, the initial impulsive phase of the EK Dra flare was not captured, and the 
key diagnostic of the hot-up/cool-down scenario - [Fe XXI] - was recorded only during the 
later phases of the event decay (although note in Fig. 4b the apparent small blueshifts of the 
coronal forbidden line compared with the average visits 3-5 behavior, which would be more 
exaggerated with the Young et ah reference wavelength). 

It also must be acknowledged that the peculiarities of the EK Dra flare decay should be 
viewed with some caution, because rarely has such an event been recorded in as much detail 
as provided by the GOS part of the HST campaign: the oddities of the emission Doppler 
shifts might well be associated with an unusual spatial orientation of the event, for example. 

In the end, highly exaggerated stellar incidents like the EK Dra FUV flare dutifully 
cause us to look back to the Sun in an effort to identify associations with solar high-energy 



transient phenomena, snch as the coronal rain mentioned earlier. Nevertheless, drawing 
snch conclusions in general is tricky given the rather tepid activity levels displayed by our 
local star, even at its worst. Probably a good thing for habitability of Planet Earth, but a 
continual challenge for the solar-stellar connection. 
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